dblogr.com/

Crop Production In The Netherlands

Graphs of crop production in the Netherlands with FAO data

Prepare Data

# devtools::install_github("derekmichaelwright/agData")
library(agData) 
library(ggforce)
library(gganimate)
myCaption <- "www.dblogr.com/ or derekmichaelwright.github.io/dblogr/ | Data: FAOSTAT"

PDF - All Netherland Crops

# Prep data
myColors <- c("darkgreen", "darkred", "darkgoldenrod2")
xx <- agData_FAO_Crops %>% 
  filter(Area == "Netherlands") %>%
   mutate(Value = ifelse(Measurement %in% c("Area Harvested", "Production"),
                        Value / 1000000, Value / 1000),
         Unit = plyr::mapvalues(Unit, c("hectares", "tonnes", "kg/ha"), 
                  c("Million Hectares", "Million Tonnes", "Tonnes/ Hectare")))
myCrops <- unique(xx$Item)
# Plot
pdf("figures_crops_netherlands_fao.pdf", width = 12, height = 4)
for(i in myCrops) {
  xi <- xx %>% filter(Item == i)
  print(ggplot(xi, aes(x = Year, y = Value, color = Measurement)) +
    geom_line(size = 1.5, alpha = 0.7) +
    facet_wrap(. ~ Measurement + Unit, scales = "free_y", ncol = 3) +
    scale_color_manual(values = myColors) +
    scale_x_continuous(breaks = seq(1960, 2020, by = 5) ) +
    theme_agData(legend.position = "none", 
                 axis.text.x = element_text(angle = 45, hjust = 1)) +
    labs(title = i, y = NULL, x = NULL, caption = myCaption) )
}
dev.off()

All Netherland Crops

# Prep data
xx <- agData_FAO_Crops %>% 
  filter(Area == "Netherlands", Measurement == "Production")
# Plot
ggcropplot <- function(x) {
  ggplot(xx, aes(x = Year, y = Value / 100000)) + 
  geom_line(size = 1, color = "darkgreen", alpha = 0.7) + 
  scale_x_continuous(limits       =   c(1960, 2020),
                     breaks       = seq(1960, 2020, 20), 
                     minor_breaks = seq(1960, 2020, 10)) +
  theme_agData() +
  theme(legend.position = "none") +
  labs(caption = myCaption, y = "Tonnes (x 100,000)", x = NULL)
}
# Plot
mp1 <- ggcropplot(xx) + 
  facet_wrap_paginate(Item ~ ., scales = "free_y", labeller = label_wrap_gen(width = 45), 
                      ncol = 5, nrow = 5, page = 1) +
  labs(title = "Crop Production - Netherlands (1/2)")
mp2 <- ggcropplot(xx) +
  facet_wrap_paginate(Item ~ ., scales = "free_y", labeller = label_wrap_gen(width = 45),
                      ncol = 5, nrow = 5, page = 2) +
  labs(title = "Crop Production - Netherlands (2/2)")
ggsave("crops_netherlands_01.png", mp1, width = 15, height = 8)
ggsave("crops_netherlands_02.png", mp2, width = 15, height = 8)

Netherlands vs Europe

# Prep data
myCrops <- c("Potatoes", "Sugar beet", "Wheat", "Tomatoes", "Carrots and turnips",
             "Onions and shallots, dry (excluding dehydrated)")
xx <- agData_FAO_Crops %>% 
  filter(Area %in% c("Netherlands", "Europe"), Item %in% myCrops,
         Measurement == "Yield") %>%
  mutate(Area = factor(Area, levels = c("Netherlands", "Europe")),
         Item = plyr::mapvalues(Item, myCrops[6], "Onions"))
# Plot
mp <- ggplot(xx, aes(x = Year, y = Value / 1000, color = Area)) +
  geom_line(alpha = 0.7, size = 1.5) +
  facet_wrap(Item ~ ., scales = "free_y", ncol = 3) + 
  scale_color_manual(name = NULL, values = c("darkgreen", "darkblue")) +
  scale_x_continuous(breaks = seq(1960, 2020, 10), minor_breaks = NULL) +
  theme_agData(legend.position = "bottom") +
  labs(title = "Yields in Netherlands and Europe",
       x = NULL, y = "Tonnes / Hectare", caption = myCaption)
ggsave("crops_netherlands_03.png", mp, width = 10, height = 6)

Bar Chart Racer

Production

# Prep data
fnames <- paste(", fresh n.e.c.", " \\(excluding dehydrated)",
               " \\(Capsicum spp. and Pimenta spp.)", ", crude", sep = "|")
xx <- agData_FAO_Crops %>% 
    filter(Area == "Netherlands", Measurement == "Production",
           !Item %in% c("Watermelons")) %>% 
  group_by(Item) %>%
  mutate(Value = movingAverage(Value, n = 3)) %>%
  group_by(Year) %>%
  arrange(Year, -Value) %>%
  mutate(Rank = 1:n()) %>%
  filter(Rank < 15) %>% 
  arrange(desc(Year)) %>%
  mutate(Item = gsub(fnames, "", Item)) %>%
  mutate(Item = factor(Item, levels = unique(.$Item)),
         Value = Value / 1000000)
# Plot
mp <- ggplot(xx, aes(xmin = 0, xmax = Value, 
                     ymin = Rank - 0.45, ymax = Rank + 0.45, y = Rank, 
                     fill = Item)) + 
  geom_rect(alpha = 0.7, color = "black") + 
  scale_fill_manual(values = agData_Colors) +
  scale_x_continuous(limits = c(-4,max(xx$Value)),
                     breaks = seq(0, 10, by = 1),
                     minor_breaks = seq(0, 10, by = 50.)) +
  geom_text(aes(label = Item), col = "black", hjust = 1, x = -0.1) +
  scale_y_reverse() +
  theme_agData(legend.position = "none",
               axis.text.y = element_blank(), 
               axis.ticks = element_blank(),
               panel.grid.major.y = element_blank(),
               panel.grid.minor.y = element_blank()) + 
  labs(title = paste("Netherlands - Production -", "{frame_time}"),
       x = "Million Tonnes (3 year moving average)", 
       y = NULL, caption = myCaption) +
  transition_time(Year) 
anim_save("crops_netherlands_gif_01.gif", mp,
          nframes = 600, fps = 30, end_pause = 150, 
          width = 900, height = 600, res = 150, units = "px")

Area Harvested

# Prep data
fnames <- paste(", fresh n.e.c.", " \\(excluding dehydrated)",
                " \\(Capsicum spp. and Pimenta spp.)", " and horse beans",
                ", processed but not spun", ", Fibre Equivalent", sep = "|")
xx <- agData_FAO_Crops %>% 
    filter(Area == "Netherlands", 
           Measurement == "Area Harvested") %>% 
  group_by(Year) %>%
  arrange(Year, -Value) %>%
  mutate(Rank = 1:n()) %>%
  filter(Rank < 15) %>% 
  arrange(desc(Year)) %>%
  mutate(Item = gsub(fnames, "", Item)) %>%
  mutate(Item = factor(Item, levels = unique(.$Item)),
         Value = Value / 1000000)
# Plot
mp <- ggplot(xx, aes(xmin = 0, xmax = Value, 
                     ymin = Rank - 0.45, ymax = Rank + 0.45, y = Rank, 
                     fill = Item)) + 
  geom_rect(alpha = 0.7, color = "black") + 
  scale_fill_manual(values = agData_Colors) +
  scale_x_continuous(limits = c(-0.08,max(xx$Value)),
                     breaks = seq(0, 0.2, by = 0.05),
                     minor_breaks = seq(0, 0.2, by = 0.05)) +
  geom_text(aes(label = Item), col = "black", hjust = 1, x = -0.002) +
  scale_y_reverse() +
  theme_agData(legend.position = "none",
               axis.text.y = element_blank(), 
               axis.ticks = element_blank(),
               panel.grid.major.y = element_blank(),
               panel.grid.minor.y = element_blank()) + 
  labs(title = paste("Netherlands - Area Harvested -", "{frame_time}"),
       x = "Million Hectares", y = NULL, caption = myCaption) +
  transition_time(Year) 
anim_save("crops_netherlands_gif_02.gif", mp,
          nframes = 600, fps = 30, end_pause = 150, 
          width = 900, height = 600, res = 150, units = "px")

© Derek Michael Wright